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Abstract 



This article is the written version of a talk delivered at the Bexbach Colloquium 
of Science 2000 and starts with an introduction into quantum chaos and its rela- 
tionship to classical chaos. The Bohigas-Giannoni-Schmit conjecture is formulated 
and evaluated within random-matrix theory. Several examples of physical systems 
exhibiting quantum chaos ranging from nuclear to solid state physics are presented. 
The presentation concludes with recent research work on quantum chromodynam- 
ics and the quark-gluon plasma. In the case of a chemical potential the eigenvalue 
spectrum becomes complex and one has to deal with non-Hermitian random-matrix 
theory. 

1 Classical Chaos 

Systems of classical mechanics can be divided into two classes: integrable and non- 
integrable systems. If there exists a constant of motion -F(p, q) for the Hamiltonian 
H{p^ q) this is connected with a symmetry 



For a system with n degrees of freedom with n constants of motion Fj, z = 1, n, and 
[Fi,Fj\ = 0, a canonical transformation to action-angle variables, (p, — > {J,0), can be 
performed 
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with frequencies tUj. The new equations of motions can be integrated leading to 



Ji = constant 
6i = ujit + ipi . 

Such a system is called integrable and each trajectory lies on an n-dimensional torus. If 
this procedure is not possible, we have to deal with a non-integrable systems. When one 
adds a non-integrable Hamiltonian Hi to an integrable one Hq 

H{J,d)=Ho{J)+eHi{J,e) , (4) 

one can show that the system cannot be integrated by perturbation theory for rational 
frequency ratios 

UJi r 



Uj s 



T,s e Z . (5) 



The famous KAM theorem due to Kolmogorov, Arnold, and Moser states that tori with 
^ sufficiently irrational are stable under small perturbations. Non-integrable systems 
possess trajectories filling the phase space ergodically and are intrinsically related to 
chaotic motion. With the KAM-theorem we get insight into the behavior of a system 
which is driven by some (order) parameter from regularity to chaos [Q. 



2 Quantum Chaos 

It is a fascinating question in which manner classical chaos is refiected in quantum systems. 
Quantum mechanics is associated with the time dependent Schrodinger equation 

= th^^ip . (6) 

In the stationary case 

ip{x, t) = (f){x) exp {—iojt) (7) 
one deals with the time-independent Schrodinger equation 

/\<pix) + Vix)(pix) = E(j)(x) . (8) 

2m 

This is a linear differential equation and its solution behaves regularly in time. Due to 
the Heisenberg uncertainty principle 

ApAx > I (9) 
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Figure 1: Chaotic billiards: Bunimovicli stadium (top left), Sinai billiard (top right), 
Pascalian snail (bottom left) and annular billiard (bottom right). Taken from Ref. [0. 

the concept of trajectories is not adequate. Therefore, the question has be posed: Are 
there differences in the eigenvalue spectra of classically integrable and non-integrable 
systems? 

Billiards became a preferred playground to study both the classical and quantum case. 
Some of the most important ones are depicted in Fig. |l|. It was the arrival of computers 
with increasing power in the late seventies when diagonalization of matrices with rea- 
sonable size became possible. The behavior of the distribution of the spacings between 
neighboring eigenvalues turned out to be a decisive signature. In 1979 McDonald and 
Kaufman performed a comparison between the spectra from a classically regular and a 
classically chaotic system . As seen in Fig. ^ they observed a qualitatively different be- 
havior between the nearest-neighbor spacing distribution of the circle and the stadium. In 
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Figure 2: Nearest-neighbor spacing distributions of eigenvalues for a circle (left) and the 
Bunimovich stadium (right). Taken from Ref. 

the first case the spacings are clearly concentrated around zero while they show repelling 
character in the second case. There were several authors contributing to this discus- 
sion and we mention the papers by Casati, Valz-Gris, and Guarneri 0, by Berry [§, by 
Robnik ^ and by Seligman, Verbaarschot, and Zirnbauer |^. 

Very accurate results were obtained for the classically chaotic Sinai billiard by Bohigas, 
Giannoni, and Schmit (see Fig. |^) which led them to the important conclusion [§: Spec- 
tra of time-reversal invariant systems whose classical analogues are K systems show the 
same fluctuation properties as predicted by the Gaussian orthogonal ensemble (GOE) of 
random-matrix theory (RMT). K systems are most strongly mixing classical systems with 
a positive Kolmogorov entropy. The conjecture turned out valid also for less chaotic (er- 
godic) systems without time-reversal invariance leading to the Gaussian unitary ensemble 
(GUE). 
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Figure 3: Nearest-neighbor spacing distributions of eigenvalues for the Sinai bilhard with 
the Wigner surmise compared to the Poisson distribution. The histogram comprises about 
1000 consecutive eigenvalues. Taken from Ref. [§. 

3 Random Matrix Theory 

In lack of analytical or numerical methods to obtain the spectra of complicated Hamil- 
tonians, Wigner and Dyson analyzed ensembles of random matrices and were able to 
derive mathematical expressions. A Gaussian random matrix ensemble consists of square 
matrices with their matrix elements drawn from a Gaussian distribution 

One distinguishes between three different types depending on space-time symmetry clas- 
sified by the Dyson parameter = 1, 2, 4 0. The Gaussian orthogonal ensemble (GOE, 
jSo = ^) holds for time-reversal invariance and rotational symmetry of the Hamiltonian 

Hmn Hjim ^nm " ( ) 

When time-reversal invariance is violated and 

one obtains the Gaussian unitary ensemble (GUE, jSo = 2). The Gaussian symplectic 
ensemble (GSE, Pd = 4) is in correspondence with time-reversal invariance but broken 
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rotational symmetry of the Hamiltonian 

HjSlh - I E , (13) 

7=1 

with H^^^ real and symmetric and if*^'^-' real and antisymmetric. 

The functional form of the distribution P{s) of the neighbor spacings s between consec- 
utive eigenvalues for the Gaussian ensembles can be approximated by 

Pp^ (s) = ap^ exp [-bp^ ) , (14) 

which is known as the Wigner surmise and reads for example in the case /5o = 2 (GUE) 

p(^s) = ^s^exp(--sA . (15) 
vr^ \ n J 

If the eigenvalues of a system are completely uncorrelated one ends up with a Poisson 
distribution for their neighbor spacings 

P{s) = exp (-s) . (16) 

An interpolating function between the Poisson and the Wigner distribution is given by 
the Brody distribution reading for the GOE case 

P(s,u;) = a(cj + l)s'^exp(-as'^+i) , " = "^"^^ (^^) ' < < 1 . (17) 

4 Examples of Physical Systems 

So far the mathematical relationship between classically chaotic systems and their eigen- 
values after quantization. Remarkably, the Wigner distribution could be observed in a 
number of systems by physical experiments and computer simulations evading the whole 
quantum world 0. 



4.1 Atomic Nuclei 



An impressive manifestation of the Wigner surmise came from nuclear physics [|T0[- The 



histogram in Fig. ^ shows the distribution of spacings of nuclear levels versus the variable 
s, the actual spacing in units of the mean level spacing D. The data set comprises 1726 
spacings of levels of the same spin and parity from a number of different nuclei. These data 
were obtained from neutron time-of-flight spectroscopy and from high-resolution proton 
scattering. The similarity to Fig. 13 is striking. 
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Figure 4: Nearest-neighbor spacing distribution for a "nuclear data ensemble" (NDE) 
compared to the RMT prediction labeled GOE and the Poisson distribution. Taken from 
Ref. pi. 
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4.2 Hydrogen Atom in a Magnetic Field 

Rydberg levels of the hydrogen atom in a strong magnetic field have a spacing distribution 
which once again agrees with RMT, see Fig. ^ [0. The levels are taken from the vicinity 
of the scaled binding energy E. Solid and dashed lines are fits, except for the bottom 
figure which represents the GOE. Additionally, a transition from Poisson to GOE behavior 
with increasing field strength is clearly visible here. 



4.3 Metal-Insulator Transition 

The conventional way to describe the localization transition is to use the Anderson Hamil- 
tonian [12 

H = '^eiajai - ^a]ai , (18) 

where aj and are the electron creation and annihilation operators at site i, subscript 
j denotes adjacent to i, and is the random energy of site i in units of the overlap 
energy of adjacent sites and is uniformly distributed in the range from [—W/2, W/2]. The 
eigenstates of the Anderson Hamiltonian in the vicinity of the energy e = experience a 
localization transition with increasing W. For a simple cubic lattice, the "metal-insulator" 
transition occurs ai W = Wc = 16 ± 0.5 and exhibits a crossover from Wigner to Poisson 
behavior, see Fig. |^. We remark that the Anderson model together with the Hubbard 
model serve also as fundamental theories for high-temperature super-conductivity. 

4.4 Yang-Mills-Higgs System 

The Lagrangian density of the SU(2) Yang-Mills-Higgs (YMH) system is given by 



£™ = i(D,0) + (Z}^'0) - V{<P) - , (19) 



where the minimal coupling of the scalar field to the gauge field and the field 
strength tensor can be written as 

iD^<P) = d,<P-tgA'^T'<P, (20) 

= 5.^^ - d.A^, + g^^'^^Al , (21) 



8 




Figure 6: Nearest-neighbor spacing distribution P{s) for the 3d Anderson model at dif- 
ferent W. Wigner and Poisson functions are also shown for comparison. Taken from 



Ref. 13 



with = c7*/2, b = 1,2,3, generators of the SU(2) algebra. The potential of the scalar 
field (the Higgs field) is 

r(0) = ^2|0p + A|0r- (22) 
This work is in the (2+l)-dimensional Minkowski space (yU = 0, 1, 2) and chooses spatially 



homogeneous Yang-Mills and Higgs fields [|1J]. When /i > the potential V has a 




minimum at |0| = 0, but for /i^ < the minimum is at 

which is the non-zero Higgs vacuum. This vacuum is degenerate and after spontaneous 
symmetry breaking the physical vacuum can be chosen 0o = (0, 0, v). If A\ = qi, = q2 
and the other components of the Yang-Mills fields are zero, in the Higgs vacuum the 
Hamiltonian of the system reads 

H=\{pI + pI) + g^v\ql + ql) + \g\lql , (24) 

where pi = qi and P2 = q2- Here w"^ = 2g'^v'^ is the mass term of the Yang-Mills fields. This 
YMH Hamiltonian is a toy model for classical non-linear dynamics, exhibiting a classical 
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Figure 7: P{s) distribution for the homogeneous Yang-Mills-Higgs system at (a) v = 1 
{io = 0.92), (b) V = 1.1 {uj = 0.34) and (c) v = 1.2 {u = 0.01), where oo is the Brody 
parameter in Eq. (|l^. First 100 energy levels and interaction g = I. The dotted, dashed 
and solid curves stand for Wigner, Poisson and Brody distributions, respectively. Taken 
from Ref. 
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chaos-order transition. To outline its connection to the quantal fluctuations of the energy 
levels, we plot in Fig. |^ the P{s) distribution for different values of the parameter v. The 
figure shows a Wigner-Poisson transition by increasing the value v of the Higgs field in 
the vacuum. By using the P{s) distribution and the Brody function it is possible to give 
a quantitative measure of the degree of quantal chaoticity of the system. In this study the 
influence of the Higgs coupling on the gauge field is analyzed for a spatially homogeneous 
YMH system. In the next section we will address the dependence of the gauge coupling 
and physical temperature on the quark field in QCD. 



5 Quantum Chromodynamics 

The Lagrangian of quantum chromodynamics (QCD) consists of a gluonic part 

Cq^^ and a part £p'"° from the quarks 



£QCD ^ ^QCD _^ ^QCD 



= -iF;,(x)Fr(x) + Y: mf)^Pfi^) , (25) 

with the Dirac spinor tpj, the quark mass rrif, the number of flavors Nf, and the generalized 
field strength tensor 

Frix) = d^A:{x) - d'^A^.ix) - gU,A^,{x)A:{x) , (26) 

where the gauge field A^ with the SU(3) indices a,b,c = 1, . . . , 8, the coupling constant g 
and the structure constants fate of SU(3) enter. The main object of study is the eigenvalue 
spectrum of the Dirac operator of QCD in 4 dimensions 

p = ^ + zg^^^ = 7^5, + 2^77M^;y . (27) 

with the Aq the generators of the SU(3) color-group (Gell-Mann matrices). Discretizing 
the Dirac operator on a lattice in Euclidean space-time and applying the Kogut-Susskind 
(staggered) prescription, leads to the matrix 



{Mks)Tx' - ^ H [Sx+ii,x' ^x^l U^^ - 5x^x'+(i ^x'f, Ul^" , (28) 

where 



U^^ = exp(tgA;ix)^) (29) 
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are the gauge field variables on the lattice and r^;^ a representation of the 7^-matrices. 

In the following we report on work of our own, partly in collaboration with B.A. Berg, 
M.-P. Lombardo, and T. Wettig. In RMT, one has to distinguish several universality 
classes which are determined by the symmetries of the system. For the case of the QCD 
Dirac operator, this classification was done in Ref. |jl5|. Depending on the number of 



colors and the representation of the quarks, the Dirac operator is described by one of 
the three chiral ensembles of RMT. As far as the fluctuation properties in the bulk of 
the spectrum are concerned, the predictions of the chiral ensembles are identical to those 
of the ordinary ensembles in Sect. H [|16[. In Ref. |T^, the Dirac matrix was studied 



for color-SU(2) using both Kogut-Susskind and Wilson fermions which correspond to the 
chiral symplectic (chSE) and orthogonal (chOE) ensemble, respectively. Here [TE\, we 



additionally study SU(3) with Kogut-Susskind fermions which corresponds to the chiral 
unitary ensemble (chUE). The RMT result for the nearest-neighbor spacing distribution 
can be expressed in terms of so-called prolate spheroidal functions, see Ref. |T^. A very 



good approximation to P{s) is provided by the Wigner surmise for the unitary ensemble, 

P^{s) = ^s'e-'^'/-. (30) 



TT 



We generated gauge field configurations using the standard Wilson plaquette action for 
SU(3) with and without dynamical fermions in the Kogut-Susskind prescription. We 
have worked on a 6^ x 4 lattice with various values of the inverse gauge coupling j3 = 
both in the confinement and deconfinement phase. We typically produced 10 independent 
equihbrium configurations for each /3. Because of the spectral ergodicity property of RMT 
one can replace ensemble averages by spectral averages if one is only interested in bulk 
properties. 

The Dirac operator, p = <f) + ig^, is anti-Hermitian so that the eigenvalues A„ of ip are 
real. Because of {^,75} = the non-zero A„ occur in pairs of opposite sign. All spectra 
were checked against the analytical sum rules J2n = ^ and Z]a„>o = 3^, where V 
is the lattice volume. To construct the nearest-neighbor spacing distribution from the 



eigenvalues, one first has to "unfold" the spectra [19 



Figure ^ compares P{s) of full QCD with Nf = 3 flavors and quark mass ma = 0.05 to 
the RMT result. In the confinement as well as in the deconfinement phase we observe 
agreement with RMT up to very high (3 (not shown). The observation that P{s) is not 



influenced by the presence of dynamical quarks is expected from the results of Ref. |]l| 



which apply to the case of massless quarks. Our results, and those of Ref. |I3 , indicate 
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Confinement = 5.2 
ma = 0.05 



Deconfinement [3 = 5.4 
ma = 0.05 




Figure 8: Nearest-neighbor spacing distribution P{s) for tlie Dirac operator on a 6^ x 4 
lattice in full QCD (histograms) compared with the random matrix result (solid lines). 
There are no changes in P{s) across the deconfinement phase transition. 




12 3 



s 

Figure 9: Nearest-neighbor spacing distribution P{s) for the free Dirac operator on a 
53 X 47 X 43 X 41 lattice compared with a Poisson distribution, e~*. 
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Figure 10: Nearest-neighbor spacing distribution P{s) for U(l) gauge theory on an 8^ x 6 
lattice in the confined phase (left) and in the Coulomb phase (right). The theoretical 
curves are the chUE result, Eq. (0), and the Poisson distribution, Pp(s) = exp(— s). 

that massive dynamical quarks do not affect P{s) either. 

No signs for a transition to Poisson regularity are found. The deconfinement phase transi- 
tion does not seem to coincide with a transition in the spacing distribution. For very large 
values of f3 far into the deconfinement region, the eigenvalues start to approach the degen- 
erate eigenvalues of the free theory, given by = J2fi=i sin^(27rn^/L^)/a^, where a is the 
lattice constant, is the number of lattice sites in the //-direction, and = 0, . . . , L^ — 1. 
In this case, the nearest-neighbor spacing distribution is neither Wigner nor Poisson. It is 
possible to lift the degeneracies of the free eigenvalues using an asymmetric lattice where 
Lx, Ly, etc. are relative primes and, for large lattices, the distribution is then Poisson, 
Pp{s) = e~^, see Fig. |[ 

We have also investigated the staggered Dirac spectrum of 4d U(l) gauge theory which 
corresponds to the chUE of RMT but had not been studied before in this context. At 
(3c ~ 1-01 U(l) gauge theory undergoes a phase transition between a confinement phase 
with mass gap and monopole excitations for (3 < (3c and the Coulomb phase which exhibits 
a massless photon for (3 > (3c [§0|- for SU(2) and SU(3) gauge groups, we expect the 
confined phase to be described by RMT, whereas free fermions are known to yield the 
Poisson distribution (see Fig. |]). The question arose whether the Coulomb phase would 
be described by RMT or by the Poisson distribution The nearest-neighbor spacing 
distributions for an 8^ x 6 lattice at /3 = 0.9 (confined phase) and at /3 = 1.1 (Coulomb 
phase), averaged over 20 independent configurations, are depicted in Fig. |10|. Both are 
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consistent with the chUE of RMT. 



Physical systems which are described by non-Hermitian operators have attracted a lot of 
attention recently, among others QCD at non-zero chemical potential fi ||2^. A formula- 
tion of the QCD Dirac operator at /i 7^ on the lattice in the staggered scheme is given 
by H 



M.,y{U, /i) 



h.c.l + 



^xt'- xt^ "y,x-H 



(31) 



u=x,y,z 



with the link variables U and the staggered phases F. For yU > 0, the Dirac operator loses 
its Hermiticity properties so that its eigenvalues become complex. The aim of the present 
analysis is to investigate whether non-Hermitian RMT is able to describe the fluctuation 
properties of the complex eigenvalues of the QCD Dirac operator. The eigenvalues are 
generated on the lattice for various values of /i. We apply a two-dimensional unfolding pro- 



cedure p4[ to separate the average eigenvalue density from the fluctuations and construct 
the nearest-neighbor spacing distribution, P{s), of adjacent eigenvalues in the complex 
plane. Adjacent eigenvalues are defined to be the pairs for which the Euclidean distance 
in the complex plane is smallest. The data are then compared to analytical predictions 
of the Ginibre ensemble |^5|] of non-Hermitian RMT, which describes the situation where 
the real and imaginary parts of the strongly correlated eigenvalues have approximately 
the same average magnitude. In the Ginibre ensemble, the average spectral density is 
already constant inside a circle and zero outside. In this case, unfolding is not necessary, 
and -P(s) is given by 



Pg{s) = cp{cs) , p{s) 



2 s lim 

Af^oo 



■7V-1 

n ^nis' 
.n=l 



N-1 

E 

n=l 



-,2n 



nler, 



(32) 



where en{x) = Ylm=o x"^/^^- and c = ds sp{s) = 1.1429.... For uncorrected eigenvalues 
in the complex plane, the Poisson distribution becomes |2^ 



s e 



(33) 



This should not be confused with the Wigner distribution (|30|) . 



We report on simulations done with gauge group SU(2) on a 6^ lattice using P = A/g'^ = 1.3 
in the confinement region for Nf = 2 flavors of staggered fermions of mass ma = 0.07. For 



this system the fermion determinant is real and lattice simulation become feasible 



We sampled 160 independent configurations |2^. In the case of color-SU(2), the staggered 
Dirac operator has an extra anti-unitary symmetry [p9| and falls in the symmetry class 
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ma = 0.07 
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/i = 0.4 
ma = 0.07 




Figure 11: Nearest-neighbor spacing distribution in the complex plane for two-color QCD 
with /i in the confinement (left) and deconfinement (right) phase. The short-dashed curve 
is the Wigner distribution for the chSE and the solid curve is the Ginibre distribution of 
Eq. m- 



with Dyson parameter /3d = 4: |^. However, one can show that the nearest-neighbor 
spacing distribution in the bulk of the spectrum for this class is also given by Eq. (P3). 



Our results for P{s) are presented in Fig. |TT[ As a function of the chemical potential /x, 
we expect to find a transition from Wigner to Ginibre behavior in P{s). This was clearly 
seen in color-SU(3) with A'^^ = 3 fiavors and quenched chemical potential p4|, where 
differences between both curves are more pronounced. For the symplectic ensemble of 
color-SU(2) with staggered fermions, the Wigner and Ginibre distributions are very close 
to each other and thus harder to distinguish. They are reproduced by our preliminary 
data for /X = and /i = 0.4, respectively. 



6 Conclusion 

We have outlined the universal applicability of random-matrix theory. Several examples of 
quantum chaos from the literature have been shown , both from numerical simulation 
and physical experiment. In some of these examples the transition from regularity to 
chaoticity could be observed, both for the classical and quantum system. 
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Concerning our own studies of quantum chromodynamics, we were able to demonstrate 
that the nearest-neighbor spacing distribution P{s) of the eigenvalues of the Dirac oper- 
ator agrees perfectly with the RMT prediction both in the confinement and quark-gluon 
plasma-phase. This means that QCD is governed by quantum chaos in both phases. We 
could show that the eigenvalues of the free Dirac theory yield a Poisson distribution related 
to regular behavior. Our investigations tell us that the critical point of the spontaneous 
breaking of chiral symmetry does not coincide with a chaos-to-order transition. 
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